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Abstract. It is well-known that Bowen-York initial data contain spurious radiation. 
Although this "junk" radiation has been seen to be small for non-spinning black-hole 
binaries in circular orbit, its magnitude increases when the black holes are given spin. 
It is possible to reduce the spurious radiation by applying the puncture approach to 
multiple Kerr black holes, as we demonstrate for examples of head-on collisions of 
equal-mass black- hole binaries. 
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1. Introduction 

Fully general relativistic numerical simulations of black-hole binaries have one 
remarkable advantage over all analytic approximation techniques (perturbation theory, 
close-limit analyses and post-Newtonian calculations): given the correct initial data, 
simulations take into account all of the physics in the problem. One may be limited by 
computer resources and the accuracy of a given numerical technique, but in principle 
no physics has been oversimplified or overlooked in the calculation. 

We are left with one qualifier: the initial data must be physically correct. Producing 
good initial data involves finding physically realistic solutions of the constraint equations 
of general relativity, and providing the appropriate physical parameters for those 
solutions. The most commonly used initial data for recent evolutions of black-hole 
binaries have been puncture data based on the Bowen- York extrinsic curvature [TJ [2] 
- of the eight codes reported to have successfully performed long-term evolutions of 
several orbits [3j HJ El El [7J El El HQ] , six use Bowen- York puncture data. 

One problem with these data is that they contain unwanted spurious gravitational 
radiation. For example, a spinning Bowen- York black hole is a Kerr black hole plus 
a Brill wave [TT]. In evolutions of black- hole binaries the spurious radiation can be 
recognized as an initial "bump" in the extracted gravitational waveform. It is interesting 
that this bump also appears in evolutions performed using data produced within the 
conformal thin sandwich (CTS) formulation of the initial- value equations |12j . and is 
even present in evolutions that start with boosted scalar fields that later collapse to 
form black holes [3]. 

In Bowen- York data, we know that the initial burst of spurious radiation is partly 
due to the Bowen- York form of the conformal extrinsic curvature. The magnitude of 
the spurious radiation grows with both the linear momentum and spin assigned to each 
black hole. We may therefore worry that the "junk" radiation will significantly perturb 
the desired orbit of the binary. For unequal-mass binaries this worry is justified: one 
sees that the spurious radiation imparts a significant initial "kick" to the black holes 

In this paper we describe a technique to reduce this spurious radiation for spinning 
black holes, using a superposition of Kerr solutions in puncture form. This technique 
is similar in spirit to the superposition of Kerr-Schild black holes [TH [To] [T6l [T7] . but 
with the advantage that it easily fits into the puncture framework that has proved so 
useful for recent black-hole simulations. In the following sections we will summarize the 
construction of black-hole initial data, outline the new method to produce superposed 
Kerr punctures, and present first results of head-on collisions of Kerr punctures, which 
demonstrate that the junk radiation has indeed been reduced. 
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2. Decomposition of the initial-value equations 

We start with the standard 3+1 decomposition [18], whereby the Einstein equations 
are separated into initial-value and evolution equations. An additional conformal 
decomposition results in initial-value equations for which there are free data, which 
must be specified, and constrained data that must satisfy the conformally decomposed 
constraint equations. 

The Bowen- York solution to the momentum constraint was originally found in the 
context of the conformal transverse-traceless (CTT) decomposition [Tj. This solution 
has the great advantage that the linear momenta and spins of multiple black holes can be 
easily specified; the global properties of the data are free parameters. On the other hand, 
the choices of the free data in the CTT decomposition necessary to obtain the Bowen- 
York solution were made for convenience, and without strong physical motivation. It 
has been argued [19] that better physically motivated choices can be made for the free 
data in the more recent conformal thin-sandwich decomposition [20], and that these 
may lead to more astrophysically realistic initial data. 

However, as we will see below, the most obvious physical choices of the free data in 
the CTS decomposition are the same as those used to construct the Bowen- York solution, 
and that the difference between the two constructions may largely be a difference in 
coordinate choice. 

In the 3+1 decomposition of Einstein's equations, the constraints are coupled 
elliptic equations for the spatial metric 7^ and extrinsic curvature Kij on one time 
slice: 



where Vj is the covariant derivative compatible with the spatial metric 7^. The 
constraint equations can be further decomposed using a conformal transformation. The 
spatial metric 7^ is related to a background metric 7^ via a conformal factor ip, 



The extrinsic curvature is also decomposed. In the conformal thin-sandwich 
decomposition [20], we define Uy = d t jij. The extrinsic curvature is split into its trace 
and tracefree parts, and can be written in terms of the lapse function N and the shift 
vector j3\ 



Vj(K ij -j ij K) =0, 
R + K 2 - = 0, 



(1) 
(2) 



(3) 



= tp- 2 Aij + ip%jK 



(4) 
(5) 




The conformal weights of the variables in these equations are 



(6) 



(7) 
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K — K, 



F = F, 



(8) 
(9) 
(10) 



and we are free to choose the function N. 

Let us now make the standard "quasi-equilibrium" choices of the free data [19] that 
have been used when numerically producing CTS data [211 [22l [19l [231 Ell ESI EZl EHl 
[29], [30]: conformal flatness, 7^ = /y, maximal slicing, K = 0, and initial stationarity of 
the conformal metric, = 0. If, in addition, we choose that N = 1/2, the Hamiltonian 
and momentum constraints reduce to the same form as in the conformal transverse- 
traceless decomposition: 



The Bowen- York extrinsic curvature is a solution to this form of the momentum 
constraint, and we see that Bowen- York data satisfies all of the quasi-equilibrium 
requirements of CTS data. 

The one extra ingredient in the standard construction of CTS data is another 
quasi-equilibrium requirement, dtK = 0, which yields an elliptic equation for N. This 
"extended CTS" system can be shown to admit non-unique solutions in some cases 
[3TI [32l [33] , but these problems have not prevented the construction of initial data for 
black-hole binaries. The extra requirement dtK = will make some change to the 
physics of the system, but it is not clear what that change will be, and it is certainly 
unclear whether the resulting data will be any more "astrophysically realistic" than 
Bowen- York data. Initial-data studies that calculate the innermost stable circular orbit 
(ISCO) for standard CTS data [221 ESI EH [29] find that its location better matches post- 
Newtonian predictions than for Bowen- York data [Ml [35]. However, recent numerical 
simulations of inspiralling black holes have cast doubt on how much sense it makes to 
discuss an ISCO for comparable-mass black holes, so it is not clear how much weight 
to give to the conclusions from these initial-data studies. One thing is certainly clear 
from dynamic numerical simulations: moving from Bowen- York to CTS data does not 
noticeably improve the spurious radiation content of the initial data. (Compare, for 
example, Figure 2 in [36] with Figure 2 in |12j.) 

In the following, we will choose a function N that is convenient for the numerical 
solution of the constraints. The solutions of the conformal thin-sandwich initial-value 
equations, together with the choice of N, provide us with a lapse N = ip 6 N and a shift 
vector but these are not necessarily ideal for dynamical evolutions. In particular, 
our choice for N will imply a lapse function that blows up at the punctures, and would 
produce slices that hit the singularities in an infinitely short time! For evolutions we 
instead make the standard choices N = 1 and f3 l = on the initial slice. 



VjA v = 
V 2 ^ = - 



(11) 
(12) 
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3. Puncture approach 

Given the Bowen- York solution of the momentum constraint, we are left to solve the 
Hamiltonian constraint to complete the specification of the initial data. It is convenient 
to choose the ansatz [2], 



where the m,i parametrize the masses of the black holes, and the r, are the coordinate 
distances to each "puncture". The Hamiltonian constraint now reduces to an elliptic 
equation for the function u. In this work, the Hamiltonian constraint is solved using 
the fixed-mesh-refinement multigrid solver in the BAM code [9]. 

There are many advantages to this approach. The n black holes are constructed 
with a topology that consists of n+1 asymptotically flat ends connected by n wormholes, 
and the slice, by construction, avoids the singularities. As a result, the initial-value 
equations can be solved on all of -R 3 , without any need to excise a region around each 
black hole. A second advantage is that it is possible to stably evolve these data, and 
although connection with the extra asymptotically flat ends may be lost with particular 
gauge choices [37], the slices remain singularity avoiding, and one need never resort to 
excision techniques. 

When constructing Kerr-like punctures, we will also have to numerically solve the 
momentum constraint. Although the conformal extrinsic curvature will diverge at the 
punctures, we can ensure that the elliptic solver deals only with non-divergent functions 
by a suitable choice of the conformal lapse, N. 

4. Head-on collisions with spinning punctures 

Let us first consider the spurious radiation content of spinning Bowen- York punctures. 
As a simple illustration, we will present evolutions of head-on collisions of equal-mass 
black holes with their spins directed away from each other. In this way, we expect the 
final black hole to be non-spinning, i.e., a Schwarzschild black hole. 

Evolutions are performed with the BAM code [9], using the BSSN formulation of 
the evolution equations [381 [39]. The lapse function and shift vector are initially N = 1 
and j3 l = 0, and are evolved using the 1+log [40] and T-driver [41] slicing conditions. 

Figure Q] shows the results from runs that started with punctures of masses 
Mi = M 2 = 0.5 placed at z = ±4M (where M = M 1 + M 2 ), and given spins 
jiS^I = 0, 0.6, 0.8, 0.9M 2 , such that the spins were directed away from the origin. The 
figure shows the (/ = 2,m = 0) mode of the real part of r\l/ 4 , extracted at r = 20M. 
(See [9] for a description of the wave-extraction algorithm in the BAM code.) 

The gravitational radiation from the merger is clearly visible after about 60M of 
evolution. Prior to that, some spurious radiation can be seen even when S z = 0. This 
is presumably due to some unphysical component to the Brill-Lindquist solution when 
the black holes are close. We do not expect the Brill-Lindquist solution to be physically 




(13) 



Beyond Bowen- York initial data 



6 



realistic unless the black holes are infinitely far apart, and indeed the spurious radiation 
decreases as the black holes are placed further apart, as is clearly seen in [TJ. (See also 
[4"2] for a post- Newtonian extension of Brill-Lindquist data.) 

As spin is added, extra spurious initial radiation is seen. The extra spurious 
radiation has roughly the quasinormal mode frequency of each of the individual black 
holes. For the case where \S Z \ = 0.9M 2 , the initial "junk" radiation is of comparable 
magnitude to the radiation from the black holes' merger. (We should point out that 
this is an exaggeration over the more realistic case of the merger of orbiting binaries, 
for which the merger waveform is far stronger; see |43j.) 
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Figure 1. The spurious gravitational radiation from the head-on collision of equal- 
mass (Mi = AI-2 = 0.5) Bowen- York punctures. The punctures are placed at z = ±4M 
and given spins a/Mi = \S Z \/Mf = 0,0.6,0.8,0.9, such that the spins are directed 
away from the origin. 



5. Kerr-like data for black-hole binaries 

Now we wish to reduce the initial radiation content for spinning punctures. Motivated 
by Dain's work in [HJ H5] , we extend work done by Krivan and Price in axisymmetry 
|46j . We start by writing the Kerr metric for a black hole with mass m and angular 
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momentum a /mm. quasi-isotropic coordinates [TT| I 

Jij = Sy + a 2 h(m, a, r, 9) Vy, (14) 
{(y 2 ,-xy,0),(-xy,x 2 ,0), (0,0,0)}, (15) 



4- = If + If + -~(L/^, (19) 



v 'm^ — a' 

V> = H h MKerr (16) 

2r 

4 = ^ err - (17) 

The analytic functional forms of A^ e " and wkcit are given in [TT] and [16]. 

This form of the metric, which appears in [H] , can be generalized to two black 
holes by superposition: 

7ij = hi + alhvlj + a^v^, (18) 

1 

2iV 

The background metric is a piece of the free data, so a straightforward superposition is 
acceptable. The conformal extrinsic curvature must satisfy the momentum constraint, 
so we consider it as a superposition of the conformal extrinsic curvatures from two single 
Kerr punctures, If and If, plus a correction that we must solve for numerically, and 
which we write in a form motivated by (Q. We continue to make the choice of maximal 
slicing, K = 0. This construction is similar in spirit to that applied to Kerr-Schild data 
by Matzner, Huq and Shoemaker [T3]. Another related work is [IT], where the authors 
adopted the Kerr form of the extrinsic curvature, but retained conformal flatness. 
The constraints now take the form 

V 2 ^ = - V 7 4^' + \tI>R, (20) 

o o 

A L /3* - (l(3)^V j \nN = - 2#V i (i£ 1 + A% 2 ). (21) 

Note that in the Kerr case the conformal Ricci scalar R is non-zero and diverges at the 
punctures. However, the Laplacian of ip with respect to the background metric (fTBl) 
diverges in the same way as ipR, and the final right-hand side of the Hamiltonian 
constraint is finite and sufficiently well-behaved that the code is able to generate 
numerical solutions. 

To solve the Hamiltonian constraint, we write the conformal factor in the standard 
puncture form, 



^ = l + j2^r^ + u ^^ + u, (22) 

i 

where \J m 2 — a? is motivated by the form of ip for the Kerr solution in quasi-isotropic 
coordinates, Eqn. f|T6l) . 

In order to solve the momentum constraint numerically, we make a new choice for 
the conformal lapse, 

N = ^ (23) 
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The reason for this choice is that the Afj 1,2 diverge at the punctures as r 3 , but 
introducing a factor of ip~ 3 ~ r 3 makes the right-hand side of (|2ip finite at the punctures. 

The constraints were solved using the multigrid solver in the BAM code [9], and 
the initial data (■y i j,K i j) reconstructed using Q - ({TO]) . 

Figure [2] shows the results of evolving configurations equivalent to those considered 
in Figure [1] for head-on collisions of spinning Bowen- York punctures. We see that the 
spurious radiation from the initial data deviates much less from that in the non-spinning 
Brill-Lindquist case. There are some high-frequency components to the junk radiation, 
but these may be due to the low tolerance used in the elliptic solver; this requires more 
analytic and numerical study. The main result is clear: the superposed Kerr punctures 
contain dramatically less spurious radiation than their Bowen- York counterparts. 




Figure 2. Spurious radiation from head-on collisions of spinning Bowen- York and 
Kerr-like punctures, for spins of a/M = 0.6 (top-left), 0.8 (top-right), 0.9 (lower-left), 
and all three (lower-right). The quantity shown in each plot is the (I ~ 2, m = 0) 
mode of extracted at 20 M. In each plot the dotted black line corresponds to the 
radiation from a collision of nonspinning (Brill-Lindquist) punctures. In the top-left, 
top-right and lower- left plots, the dashed line shows the radiation from Bowen- York 
punctures, as also shown in Figure [TJ and the solid line shows the radiation from Kerr- 
like punctures. The lower-right plot allows a comparison between the three choices of 
spin. 
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6. Conclusions 

We have presented a method to construct initial data for two Kerr-like black holes in 
the puncture framework. This construction is based on the form of the Kerr metric 
in quasi-isotropic coordinates [TT], and following the suggestions outlined in [44], [45] 
extends the work in axisymmetry of [4"6] . 

For the simple test case of a head-on collision of two black holes, we have 
demonstrated that the new initial data contain significantly less spurious gravitational 
radiation than their popular Bowen- York counterparts. 

The results here serve as an illustration of the potential of this method. In 
future work we will explore its effectiveness in more general configurations, and extend 
the construction to spinning black holes that also possess linear momentum. Such a 
construction would be necessary to make these data a feasible alternative to Bowen- 
York data for simulations of black-hole binary inspiral. 
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